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PROCEEDINGS 

SCIENCE 


m 

We study a possible relationships between seasonal distributions of the visually observed cloud¬ 
less days (CD) and cloudless nights (CN) at Abastumani Astrophysical Observatory (41.75N, 
42.82E; Georgia) in 1957-1993. The annual variations of monthly numbers of CD and CN have 
been observed, with maximum in August for CD and in September for CN. During geomagnetic 
disturbances it is also observed the growth of number of CD in September and March (equinoctial 
months), and for CN, together with September, in June, April and February. We assume that this 
phenomenon indicates an influence of cosmic factors on cloudiness, as well as the existence of 
semiannual and possibly shorter-periodicity variations. This cosmic factor can be the manifesta¬ 
tion of different rates of the galactic cosmic rays (OCRs) flux variations in CD and CN periods. 
The influence of GCR flux on ionization of lower atmosphere and variations of density of cloud 
condensation nuclei also can be connected to the annual and seasonal changes of temperature at 
Earth surface of this region. To comprehend behaviors of the annual and semi-annual variations 
of the GCR intensity and their possible relationships with the seasonal distributions of CD and 
CN we compose and numerically solve two dimensional (2-D) time dependent transport equation 
including all important processes in the heliosphere. An analysis of experimentally observed and 
theoretically obtained results have been carried out. 


The 34th International Cosmic Ray Conference, 
30 July- 6 August, 2015 
The Hague, The Netherlands 


* Speaker. 

^A footnote may follow. 


0 Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 







Annual Variations of Cosmic Rays 


M. V. Alania 


1. INTRODUCTION 

A relationship of changes of the cloudless days (CD) and cloudless nights (CN) numbers with 
the GCR intensity variations is studied at first time based on the unique observational data in the 
Abastumani Astrophysical Observatory (Georgia) during 1957-1993. So, in the beginning we con¬ 
sider the behavior of the average GCR intensity during the yearly period based on neutron monitors 
data. Then we compare results with similar changes in CD and CN for the same period. The inter¬ 
annual changes of atmospheric processes, including cloud covering in the lower atmosphere, can 
possibly be related to the seasonal changes of absorption of solar electromagnetic radiation energy 
by the Earth’s surface, which depend on geometry of solar-terrestrial position in the heliosphere. 
During Earth’s rotation around the Sun, the interplanetary and geomagnetic field geomefry and, 
as a resulf, the influence of fhe solar corpuscular radiation on fhe evenfs on the magnetosphere 
i] and also the magnetosphere-ionosphere-atmosphere coupling change |Q|. To generate geomag¬ 
netic disturbances by solar wind, the interplanetary magnetic field (IME) and geomagnetic field 
configurafion are effecfive af equinox monfhs (March/April and Sepfember/Ocfober) [|I]], where fhe 
z componenf of fhe IME is direcfed soufhward. The high number of geomagnefic disfurbances af 
equinoxes and a small number around solstice monfhs gives semi-annual variabilify of frequency 
of appearance of fhe phenomena characferisfic for magnefosphere-ionosphere-afmosphere coupling 
processes [||]. Solar wind disfurbances cause changes nof only in fhe geomagnefic field buf modu- 
lafe GCR flux, as well. These in furn can cause variafions of ions produced by GCR and hence fhe 
densify of cloud condensafion nuclei (CCN). Thus, fhere is a possible coupling befween changes 
in GCR flux and cloud formafion processes [Q]-[Q]. High mean day-nighf femperafure in the lower 
atmosphere and the Earth’s surface can be favorable for cloudless days and nights. When a differ¬ 
ence between day and night-time temperature is comparatively large, then cloud formation condi¬ 
tions are favorable, including the influence of cosmic facfors, like GCR. Seasonal peculiarifies of 
femperafure variafions from day fo the night as well as with height for given region of the lower 
atmosphere also can change the annual and semi-annual effect of production of CCN by GCR. This 
in turn may result in different behavior of cloud formation during the day and the night. 

Increase of the CCN density produced by GCRs flux enhancemenf, where fhe wafer vapor is 
near safurafion, should sfimulafe fhe cloud formafion processes and grow fhe cloudiness [^]. De¬ 
crease of GCRs flux in similar conditions should be favorable for cloudless days and nighfs. In fhe 
lower afmosphere where fhe clouds mosfly are formed, fhe femperafure and humidify variafions are 
differenf for day and nighf. Because of fhis, fhe processes initiating cloud formafion by GCR may 
be differenf during day and fhe nighf. Togefher wifh seasonal as well as day-nighf variafions of 
afmospheric conditions, if is possible fhaf fhe impacf of cosmic facfors on cloud covering change 
by day-nighf and seasons. Togefher wifh fhe long-ferm variafions of cloud cover, which can cor- 
relafe wifh solar and geomagnefic acfivifies, as well as galacfic cosmic ray flux [§]-[§], it should 
be important to study their inter-annual/seasonal variations, which is the main purpose of this pa¬ 
per. We will examine inter-annual distribution of CD and CN for various level of geomagnetic 
activity and variations in it caused by possible cosmic factor. We also consider the inter-annual 
variations of the main cosmic factors: GCRs flux, solar wind velocify, and IME af fhe solar acfivify 
minimum (1975-1978) and maximum (1990-1991) phases. The infer-annual disfribufion of GCRs 
flux will be considered fheorefically as well. Numerically will be solved fwo-dimensional (2-D) 
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Figure 1: Inter-annual distribution of total monthly numbers of cloudless days in Abastumani Astrophysical 
Observatory in 1957-1993 at different geomagnetic disturbances; for Ap > 7 (blue circles), Ap > 12 (red), 
Ap > 20 (green) and Ap > 30 (violet). Yellow line corresponds to all cloudless days. 
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Figure 2: Inter-annual distribution of total monthly numbers of cloudless nights in Abastumani Astrophys¬ 
ical Observatory in 1957-1993 at different geomagnetic disturbances: for Ap < 12 (blue dot), Ap < 20 (red), 
Ap < 30 (green) and Ap < 40 (violet). Yellow line corresponds to all cloudless nights. 

time-dependent transport equation including all relevant processes in the heliosphere. An analysis 
of experimentally observed and theoretically obtained results have been carried out. 

2. Experimental data 

Figure [I] shows that the biggest number of cloudless days is in August. During magnetically 
relatively disturbed conditions (Ap > 7) it shifts to September, where is observed a maximal num¬ 
ber of magnetically disturbed day-nights [|^. In the Fig. |i| the increase of number of magnetically 
disturbed cloudless days is observed, both in September and March, pointing to semi-annual vari¬ 
ations modulating the annual variations of cloudless days, which indicates the influence of cosmic 
factors on cloud covering. This phenomenon is enforced by the fact that maximal number of cloud¬ 
less nights is in September (Fig. ^ and at geomagnetically comparatively quiet periods (Ap < 40) 
it is in August again, like for cloudless days. We note that in this region the maximal day-night 
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Figure 3: Superimposed over years 1975-1978 monthly changes of the solar wind velocity V and the 
interplanetary magnetic field strength B and corresponding variations in the GCR intensity based on the 
Oulu neutron monitor data. To fit one scale, the solar wind velocity was scaled down by 400 km/s and IMF 
strength by 3.5 nT. 

temperature at the surface is observed mostly in August and, as was mentioned, the number of 
cloudless day-nights is more expectable. During magnetically disturbed conditions, such annual 
distribution is modulated by semi-annual variations that can be caused by the influence of cosmic 
factors on cloud covering. We assume that one of these cosmic factors, influencing cloud covering 
processes, can be considered the inter-annual variations of GCRs flux at cloudless days and nights. 

3. Model of the annual GCR intensity variation 

We assume that as one of the conventional candidates to affect the character of cloud formation 
and cloud coverage (cloudiness and cloudless processes) can be considered GCRs. An energy 
contribution of cosmic rays in the Earth atmosphere is smaller than the energy from a total solar 
radiation. Nevertheless, cosmic ray particles are only a source of ionization in the lower atmosphere 
[^], being further a reason for cloud formation. Unfortunately, there is not yet well-known how the 
ionization produced by cosmic ray particles takes place in microphysics of the cloud formation. 

Distributions of CD and CN have clearly expressed annual variation (Fig. |T] and Fig. |)- If 
we ascribe these variations to the cloud formation owing to the ionization by cosmic ray particles 
(of energy 5-35 GeV, to which neutron monitor (NM) respond), there is possible to estimate an 
expected annual changes of cosmic ray flux solving Parker’s transport equation of anisotropic dif¬ 
fusion 1^. We believe that any model describing annual changes of the GCR intensity giving an 
opportunity to estimate a coupling between cosmic ray flux and the cloud formation is extremely 
useful. In [^] we have solved Parker’s transport equation installing in equation the annual and semi¬ 
annual alternation of solar wind velocity comparable with experimental data. We demonstrated that 
changes in solar wind velocity in the range ~ 350-500km/s cause about 0.3% changes in the GCR 
intensity. This variation could not provide the enough changes of the level of ionization (aerosols) 
responsible for cloudiness-cloudless state of the Earth’s atmosphere [||] and references therein. 

In this paper, we compose a new 2-D time-dependent model of GCR propagation in the he¬ 
liosphere including two crucial parameters i.e. in-situ measurements of the solar wind velocity V 
and strength B of the IMF. Besides, in the present paper we consider one minimum and one max- 
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Figure 4: Superimposed over years 1990-1991 monthly changes of the solar wind velocity V and the 
interplanetary magnetic field strength B and corresponding variations in the GCR intensity based on the 
Oulu neutron monitor data. To fit one scale the solar wind velocity was scaled down by 400 km/s and IMF 
strength by 5 nT. 

imum epoch of solar activity, separately. For this purpose we performed the superposition of the 
monthly changes of the solar wind velocity V and the IMF strength B during three years in the min¬ 
imum 1975-1978 of solar activity and two years 1990-1991 in the maximum of solar activity. Fig.|| 
presents superimposed over years 1975-1978 monthly changes of the solar wind velocity V and the 
IMF strength B and corresponding variations in the GCR intensity based on the Oulu NM data. The 
same for the period of maximum of solar activity 1990-1991 is presented in Fig. The selection 
of Oulu NM with the effective rigidity of 10-12 GV is caused by: (1) Oulu NM operates stable for 
the extended period, (2) data of Oulu NM does not undergone changes due to disturbances of the 
Earth magnetosphere, (3) changes in amplitudes of various classes of the GCR intensity variations 
are relatively significant, and (4) data of Oulu NM are sensitive to any kind of global changes in 
heliosphere. 

Fig. I shows that in the minimum epoch 1975-1978 an average annual alternations (1) of solar 
wind velocity is ~ 100 km/s ( from ~ 400 up to ~ 500 km/s), (2) of the IMF strength B is ~ 1.2 
nT ( from ~ 5.05 up to ~ 6.25nT), and (3) in cosmic ray variations a range of changes is, I ~ 
1.25%. We underline that there hardly can be recognized any definite relationship (inter-annual or 
annual) among the changes of the parameters - V, B and I. Fig. ^ shows that in maximum epoch 
1990-1991 an average annual alterations (1) of solar wind velocity is ~ 90 km/s (from ~ 410 up to 
~480 km/s), (2) of the IMF strength B is ~1.0 nT ( from ~ 5.1 up to ~6.1nT), and (3) in cosmic 
ray variations a range of changes is, I ~ 11%. We underline that there hardly can be recognized 
any precise regular (inter-annual or annual) relationships among the changes of these parameters - 
V, B and I. Presented in the Fig. ^ and Fig. ^ monthly changes of the solar wind velocity V and 
IMF strength B were implemented into the Parker’s time-dependent 2-D transport equation |^] : 

^ = V.{Kfj-VN)-iv, + U)-VN+^^VU, (3.1) 

where, N and R are omnidirectional distribution function and rigidity of GCR particles, respec¬ 
tively; T-time, F-solar wind velocity, Vd is the drift velocity. We set up the dimensionless den¬ 
sity / = ^ ^ ^ = :^ time f = -^ distance; where, Nq is density in the Local Interstellar 
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Figure 5: The GCR intensity variation for the Oulu neutron monitor (red line) and expected from modeling 
for rigidity R=\A and 10 GV (magenta and blue line) results for the period of solar activity minimum. 



Figure 6: The GCR intensity variation for the Oulu NM (red line) and results expected from modeling for 
rigidity R= 14 and 10 GV (magenta and blue line) for the period of solar activity maximum. 


Medium (LISM) accepted as Nq = 4;r/o, the intensity /q in the LISM has the form: /q = 
+ 5.85r^*-^^ + T is kinetic energy in GeV; p and po are the radial dis¬ 

tance and size of the modulation region; Tq -is the characteristic time corresponding to the changes 
in heliosphere for the particular class of GCR variation. For the annual variation, we accept that Tq 
is equal to one year. A size of the modulation region is po=100AU, the upwind-downwind asym¬ 
metry of the heliosphere is not taken into account, as far po is significantly greater (more than 10 
times) than the Larmor radius of GCR particles to which NMs respond. 

The anisotropic diffusion tensor of GCR Kn = K-P consists of the symmetric K-P and 

‘7 U U 

(A) 

-antisymmetric parts. We implement a drift velocity of GCR particles as, [p^]. 

The heliospheric magnetic field vector B is faken, as B = (1 — 2H{d — B )){Br ■ + ^ 9 ^ 9 ) 0 - 

where H is fhe Heaviside sfep function changing fhe sign of fhe global magnetic field in each 
hemisphere and B' corresponds fo fhe heliolafifudinal posifion of fhe heliospheric neufral sheef 
(HNS), Cr and ^9 are fhe unite vectors directed along fhe componenfs Br and B 9 of fhe IMF for 
fhe fwo-dimensional Parker field [|I^. Parker’s spiral heliospheric magnefic field is implemenfed 
fhrough fhe angle t/r = arctan{—Bip/Br) in anisofropic diffusion tensor of GCR particles (t/r is fhe 
angle befween magnefic field lines and radial direction in fhe equatorial plane) and rafios j3 = 
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K±/K\\ and j8i = of the perpendicular K± and drift diffusion coefficients to the parallel 

^11 diffusion coefficient. A parallel diffusion coefficient used in modeling is expressed, as: A'n = 
K()K{r)K{R,a) where Kq = 1.9- K{r) = 1 + 0.5r/ro; K{R,a) = contributes to the 

changes of the parallel diffusion coefficient A'|| due to dependence on the GCR particles rigidity R. 
We assume that a = 0.7 in the minimum and a = 1.2 in the maximum of solar activity. Drift effect 
due to gradient and curvature of the regular IMF is implemented in the model by means of the ratio 
of the drift Kd diffusion coefficients to the parallel diffusion coefficient j3i = Kd/K\\. In this 
model, we consider that a drift effect during the maximum solar activity is scaled down by 30% 
(almost diffusion dominated case) with respect to the minimum of solar activity (drift dominated 
case). The equation 3.1 was transformed to the algebraic system of equations using the implicit 
finite difference scheme, and then solved by the Gauss-Seidel iteration method. The solutions for 
each layer of rigidity R (R = 100,90,80,..., 10 GV) for the stationary case are considered as an 
initial conditions for the nonstationary case for the given rigidity R and at time t = 0. The equation 


3.1| in spherical coordinate system for dimensionless variables is derived in detail in []14|], while the 


details of its numerical solution of the 3D nonstationary equation are given in []_5]. 

Results of the numerical solution of the equation ^ with included changes of the solar wind 
velocity V and IMF strength B (Fig. ^ and Fig. ^) are presented in Fig. ^ for the minimum of solar 
activity, and in Fig. ^ - for the maximum of solar activity. The expected from the model variations 
for paiticles with rigidity /? = 14 and 10 GV are compared with the changes in the GCR intensity 
observed in Oulu NM data. To obtain more realistic results we implemented into our model in 
situ measurements of V and B. It is seen from Fig. ^ and Fig. ^ that we have not achieved a very 
good agreement between observed by Oulu NM data and modeling results neither for minimum nor 
maximum epochs. However, one can recognize that there are not significant distinctions between 
range of changes for experimental data and outcomes of theoretical calculations. Of course, using 
some physically justified assumptions, e.g., the changes in parallel and perpendicular diffusion 
coefficients versus the strengths B and turbulence of the IMF, we could achieve more acceptable 
agreement between observed by Oulu NM data and modeling results. Nevertheless, our aim in this 
paper was to speculate, as less as possible and create a realistic 2-D time-dependent model based 
on in situ observations of parameters defining the general processes in cosmic ray modulation. We 
believe that CD and CN distributions are related to the processes in the cosmos. So, after careful 
study of the coupling between level of ionization in lower earth atmosphere and GCR flux, CD and 
CN could be considered as the useful proxies of the space weather conditions. 


4. Conclusions 

1. The inter-annual distribution of cloudless days and nights in Abastumani Astrophysical Ob¬ 
servatory reveals both annual and semi-annual variations for various levels of geomagnetic 
disturbances. For the annual cycle, the number of cloudless days is the biggest in August, 
which may be expectable, since in this month the mean day-night surface temperature in the 
region is maximal. In geomagnetically disturbed conditions {Ap > 7), having semi-annual 
character, the increase in a number of cloudless days at equinoxes possibly cause the shift of 
the maximal number of cloudless nights from August to September. This phenomenon and 
the fact that number of cloudless nights are the biggest in September and for less geomag- 
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netic disturbances {Ap < 40) shifts back to August, point to the influence of cosmic factors 
on cloud covering processes. 

2. One of the possible cosmic factors the OCRs flux change is considered, which reveals annual 
variations for cloudless nights, with the greatest decrease in June, where, for magnetically 
disturbed conditions, the maximal frequency of cloudless nights are observed. 

3. The OCRs flux observed by Oulu NM shows a presence of some annual and inter-annual 
variations. These inter-annual variations differ from its modeling results but gives its notice¬ 
able reduction in June and August, where the relative number of CN and the total number of 
CD are increasing, respectively. This result should be important for improving an assumed 
model (implementing into transport equation in situ measurements of V and B) for investiga¬ 
tion of the observed properties of the inter-annual variations of the CD and CN distribution. 
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